Forbidden activation levels in a non-stationary tunneling process 
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Tunneling in the presence of an opaque barrier, part of which varies in time, is investigated numerically and ana- 
lytically in one dimension. Clearly, due to the varying barrier a tunneling particle experiences spectral widening. 
However, in the case of strong perturbations, the particles' activation to certain energies is avoided. We show 
that this effect occurs only when the perturbation decays faster than t~ 2 
PACS: 73.40. Gk, 71.55 and 72.10.Fk 



Tunneling, which is a very common phenomenon in the 
quantum world, is rarely a stationary process. In most 
cases, it is involved with temporal changes of the medium 
through which the particle is tunneling. In some cases, 
such variations are mostly due to transport of the tun- 
neling particles themselves. In these cases, the problems 
are intrinsically nonlinear. However, when the source of 
the changes is external, the problem can be treated as a 
linear one. In this category, one can even find station- 
ary many-body transport problems where a single par- 
ticle can experience the influence of its surroundings (in 
this case its neighbors) as an alternating external force. 
Many works (see, for example, refs. done in this 

field have revealed an abundance of physical effects: res- 
onances and oscillations in energy, exponential increase of 
the tunneling current, activation- assisted tunneling and 
elevator resonance activation. 

This paper discusses the influence of a very strong per- 
turbation (compared with the incoming particles' energy) 
on the tunneling process. It is assumed that the per- 
turbation potential is always positive (see Fig.l), i.e., it 
increases the barrier (locally) and never creates a local 
well (inside the barrier). Therefore, effects like elevator 
resonance activation do not occur (see ref. pj). On the 
contrary, in the adiabatic approximation, the perturba- 
tion - being very strong - should reduce the transmission 
dramatically. However, when the characteristic time of 
the perturbation decreases, the particle may absorb en- 
ergy from the external perturbation, which can assist it 
in its due course of tunneling. It is clear, for example, 
that in the extreme case, where the potential barrier is 
smaller than Ti/t (t is the characteristic time of the per- 
turbation) , the particle can easily (with high probability) 
absorb enough energy to leap over the barrier. This pa- 
per discusses the intermediate case, where the process 
has a tunneling nature but allows for energy quanta to 
be absorbed by the tunneling particle. A very strong per- 
turbation, which is activated for a finite period, permits 



propagation only at the very onset and at the final decay 
of the perturbation. When the wave function in these 
two instances has the same phase, maximum transmis- 
sion is expected, but when in these two moments the 
wave function is out of phase, one should expect mini- 
mum transmission. These are the forbidden activation 
energies we are looking for. 
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FIG. 1. Activation in a tunneling process due to a vary- 
ing potential barrier. The different gray levels represent the 
temporally varying part of the potential. 

This system, which is presented in Fig.l, can be de- 
scribed in terms of the Schrodinger equation in the fol- 
lowing way 
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- -^I + U ^ + WMK*)* = i-£ (1) 

where we have used the units fi — 1 and 2m = 1 (to is 
the electron's mass), and f(x) is a function, which van- 
ishes very rapidly for x — > oo, and U{x) is the barrier's 
potential. 

We seek the solution 

iP(x,t) = ip n (x,t) + J due-^auG+ix) (2) 

where ipn(x,t) is the homogeneous solution of eq^, 
i/jfi(x,t) = ipfi(x)e~ int , while ip solves the stationary- 
state equation 



1 



d 2 ipn 
dx 2 



+ [U{x) - SI] = 0; 



(3) 



G+ is the " outgoing" Green function of equation || (but 
with a general w instead of St). 

The perturbation can be regarded as the force that is 
activated upon the particle for a finite period r (/ has 
momentum dimensions, therefore //r has dimensions of 
force) . 

A straightforward substitution of eq.^| in eq.[l] reveals 
the integral equation B 



f(u-St)pn(0) + a u 



dw'/(w-u')a W 'G ! i(0)=0 



where 



/(w) = (27T)" 1 / dtf(t/r) exp(iwt) 



(4) 



(5) 



and some tedious calculations || reveal that, for a rect- 
angular barrier with potential height ?7 and width L, 



GUO) 



coth [kL + i arctan(fc/Ac)] 
2^ 



-(2/,)" 1 (6) 



where k = y^w and ft = y/U — ui and the approxima- 
tion on the right hand side takes place for an opaque 

barrier WUL > . 

It should be noted that the results of this paper are 
valid for any opaque barrier. We have chosen a rectan- 
gular barrier since it has a relatively simple expression. 
In fact, the approximation in eqj| (which is a good es- 
timation for any opaque barrier) yields, for all practical 
purposes, the same results. 

In general, eq^ does not have an analytical solution. 
Thus, we will begin with the numerical results. In Fig. 
2 we present the results of the exact solution of eq.[| for 
the temporally confined perturbation 



/(t/r) = (/o/r)ex P (-|t/r|). 
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The x-axis represents the frequency u) (i.e., the fre- 
quency of the transmitted particles), and the y-axis 
represents the logarithm of r (the perturbation's char- 
acteristic time). Within this matrix, the gray levels 
indicate the logarithm of the transmission coefficient, 
i.e., In \a UJ G'^(x > L)\. In the figure, dark spots indicate 
a high probability of a particle being emitted with the 
corresponding frequency, or energy (u>). 

From the figure we see that, for very large r (i.e., the 
adiabatic case), the outgoing frequency is almost identi- 
cal to the incoming one, St. That is, the particles simply 
experience spectral widening. 

When t decreases (and thus the process becomes more 
energetic), the spectral width of the outgoing particles 
increases, as it should by the uncertainty principle. The 



strange part comes later, when r reaches the value To , 
after which the increase in the spectral width becomes 
more moderate. This is the point at which the acti- 
vated particles prefer to be activated to the first reso- 
nant energy ~ U + (n/2L) 2 . These resonant energies 
are barrier-dependent; therefore, one should expect dif- 
ferences in these resonances from one barrier to the next. 
However, these resonances have nothing to do with the 
following effect, which is the main issue in this paper. 

When t decreases even further, the spectral width is 
split into discrete probable frequencies-like a fan. That 
is, the transmitting (through the barrier) particles are ac- 
tivated to specific preferred frequencies, while activation 
to others frequencies is forbidden. 

To explain this behavior, we can use the approxima- 
tion that U — St r^ 1 . Since, U — SI can be arbitrar- 
ily large, this approximation still enables us to discuss 
relatively rapid changes. Within this temporal regime, 
only the spectral vicinity of SI affects the amplitude a u . 
Thus, we can approximate the Green function in eq. |] 
according to: G+(0) — * G^(0). The important thing, 
which allows for this exchange, is that when SI is very far 
from U, G+{0) ~ -{2VU^~u 



does not change much 



within the spectral range (again, recall that this approx- 
imation is almost barrier-independent). Hence, with this 
assumption in hand eq.^ can readily be solved with a 
Fourier transform: 



A{t) 



f(t/T)(p n (0)exp(-iSti) 
1-/(</t)G+(0) 



(8) 




FIG. 2. The probability of a particle being activated to 
energy u) ac t, for various characteristic times r. The y-axis 
represents r on a logarithmic scale and the x-axis represents 
the outgoing activation energy u) act . Inside the figure, the 
darker the spot the higher the probability it represents (again 
on a logarithmic scale). For an arbitrary r', the correspond- 
ing gap between adjacent forbidden activation energies, i.e., 
6u(t') is illustrated. 
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where A(t) is the inverse Fourier counterpart of a w . 
For convenience, we use the dimensionless function g: 
f(t/r) = (f /r)g(t/T) (where f = /(0)). The results in 
Fig. 2, for example, were taken for g{x) = exp(— |x|). It 
is assumed that fo = /(0) is the maximal value of the 
perturbation. Therefore, 



A(t) 



fog(t/T)cp n (0) exp(-int) 



r - f o g(t/r)G+(0) 
Eq.S determines two time scales 



Tl = /oG[ 



and 



t-2 = rg 



'H-foGa/r). 



(9) 



(10) 



(11) 



When t > Ti, the denominator of eq. ^ can be ignored, 
and the broadening of the spectral width is exactly equal 
to the spectral width of the perturbation. Note that since 
we are dealing with a strong perturbation (large /o), this 
characteristic time is very long. 

When r < Ti (and this is the more interesting regime), 
the denominator cannot be ignored, and here ti comes 
into play: when t < r 2 , A(t) ~ <pn(0) exp(-tOi)/G^(0); 
except for the oscillating term the absolute value of A{t) 
is almost a constant. However, when t 3> T2, the denomi- 
nator of eq. ^ can again be ignored (because of the rapid 
decay of g). Therefore, A{t) has a kind of rectangular 
shape: it is almost a constant for a time ~ r 2 ^> t and 
then it abruptly decays with the time scale r For such a 
"rectangular" shape, the Fourier transform is something 
like the sine function, which is oscillating and decaying 
at the same time. This oscillating "frequency" of a w is r 2 
(note that we are dealing here with the frequency space, 
hence the frequency has dimensions of time). Thus, one 
can write a u ~ a[(uj — Q,)t2\. Multiplication by the Green 
function G+(x > L) distorts the picture a little (it in- 
creases the high frequencies' contributions at the expense 
of those of the low frequencies), but it does not change 
the fact that the particles can be activated ("act") only 
by a product of an integer and the quantity Sw ~ t^ 1 . 
That is, 



O + nSu>, 



(12) 



where n is an integer. Similarly, no activation occurs 
to the forbidden ("for") activation levels 



for 



n + — ) rW . 



(13) 



These are the forbidden frequencies we were looking 
for. From the definition of r 2 (eq. |ll|), one learns that it 
depends not only on the characteristic time r and on the 
incoming frequency f2, but also on the strength of the 
perturbation /o and, maybe more surprisingly, on the 
functional behavior of /. For example, suppose 



7 (x)=exp(-|xr) 



then 



r 2 = r 



In - 



foGs 



(14) 



(15) 



It is clear that when the perturbation is strong (/o is 
large) the dependence on m can be very significant. 

Let us take the following temporal dependence that 
corresponds to a milder decay: 



9n(x) = 



1 



1 



(16) 



Such temporal behavior leads to a particularly simple 
solution, which allows tracing the onset of the " discrete 
activation level effect" . In this case 



a^ = -J dtA(t) expM) = -J d XwT -, 



(17) 



where we adopt the dimensionless parameters s = (lu - 
£1)t, a 11 = — G(0)/o/ T + 1- an d of course x = t/r. 
But eqjlT] can also be rewritten 



fo - , ^ 
-rrrffn(sa) 



(18) 



where g n is the inverse Fourier transform of g n (see 
eq.|). 

For y — » oo, an asymptotic expression can be written 



9n{y) ~ cos[?/cos(7r/n)] exp[-ysin(7r/n)], 



(19) 



(for n = 2 and n = 4 this is the only term in the 
asymptotic expansion) . 

Thus, the oscillating part of a u reads 



cos 



(uj — Q)t ( — + 1 ) CDS 



Gf 



(20) 



Similar to the exponential case (eq. |l4]) we can define 
a characteristic time 



T-2 



r(-G/ /r + l) 1/ "cos( 7 r/n) 



(21) 



which determines the activation energies according to 
eq. [l2]. Similar to eq. [l^, ti — > r when the functional be- 
havior of the perturbation dies out very quickly (n — * 00). 
Moreover, eq. ^l] also reveals the onset of the discrete ac- 
tivation pattern. The parameter n must be larger than 2, 
otherwise no oscillations occur. 

Therefore, we conclude that in order to obtain the dis- 
crete activation level effect, the perturbation must die 
out faster than (t/r)^ 2 . 

It should be noted, that because of the point-like na- 
ture of the alternating impurity, the whole derivation of 
the effect depends on the barrier via the Green function, 
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however, as approximation (g) suggests, any opaque bar- 
rier would yield similar results. 

To summarise, the activation energy in the case of a 
strongly varying potential barrier was investigated. It 
was shown that when the temporal perturbation ampli- 
tude is very large, the particles cannot be activated to 
certain discrete energies u?f or — + {m + l/2)Sui, where 
5uj ~ [t<? _1 (— /oGq/t)] -1 , t is the characteristic time 
of the perturbation, /o is its strength and g~ l is the 
inverse function of the perturbation's temporal depen- 
dence. Moreover, it was demonstrated that in order to 
see this forbidden activation effect, the perturbation must 
die out faster than i -2 . 

I would like to thank Prof. Mark Azbel for stimulating 
discussions. 
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